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Abstract. Geometric phase has found a broad spectrum of applications in 
both classical and quantum physics, such as condensed matter and quantum 
computation. In this paper we introduce an operational geometric phase for 
mixed quantum states, based on spectral weighted traces of holonomies, and we 
prove that it generalizes the standard definition of geometric phase for mixed 
states, which is based on quantum interferometry. We also introduce higher 
order geometric phases, and prove that under a fairly weak, generically satis- 
fied, requirement, there is always a well-defined geometric phase of some order. 
Our approach applies to general unitary evolutions of both nondegenerate and 
degenerate mixed states. Moreover, since we provide an explicit formula for 
the geometric phase that can be easily implemented, it is particularly well 
suited for computations in quantum physics. 



1. Introduction 

In the seminal paper pQ, Berry demonstrated that a quantum system in a pure 
state that undergoes an adiabatic cyclic evolution can pick up a phase which is 
independent of the dynamics of the system. Berry's publication began a, still ongo- 
ing, era of intensive research on holonomy effects in quantum theory, and some of 
the more important earlier papers on the topic have been reprinted in [2]. There, 
for example, one can find [3 , in which Simon identifies Berry's phase as a parallel 
transport holonomy in a Hermitian line bundle over the parameter manifold of the 
system's Hamiltonian, and [4J, in which Aharonov and Anandan extend Berry's 
result to the non-adiabatic cyclic case, and make explicit that Berry's phase is due 
to the curvature of projective Hilbert space. More recent applications of quantum 
mechanics in which holonomy effects play a crucial role include quantum informa- 
tion processing and quantum computing J5J [51 [SJ EH EH HH EH (HI • Indeed, it 
is believed that quantum computation implemented by holonomic quantum gates 
is intrinsically fault-tolerant and robust against noise [151 [TBI [T7] . 

States of experimentally prepared quantum systems generally exhibit classical 
uncertainty. Hence they are most appropriately described as mixtures of pure 
states. Uhlmann |18[ fT9] was among the first to develop a theory for geometric 
phase for parallel transported mixed states. The theory is based on the concept 
of purification. Another approach to geometric phase for mixed states, based on 
quantum interferometry, was proposed by Sjoqvist et al. [20]. They derived, and 
carefully analyzed, an expression for the intensity of the output signal in a classical 
Mach-Zehnder interferometer, and from this they extracted a definition of geometric 
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phase for mixed states. Several experiments |21[ 122 23, 24] have been carried out 
that verify Sjoqvist et al.'s results. 

A shortcoming of the definition of geometric phase given in [20] is that it only 
applies to parallel transported nondegenerate mixed states. In the current paper 
we provide a geometrical framework, inspired by arguments of Montgomery |25| . 
in which geometric phase for unitarily evolving general mixed states has a natural 
definition. We show that this definition agrees with the definition in [20] in those 
cases where their definition applies, and we derive an explicit formula for the geo- 
metric phase which can be easily implemented. Our approach is thus particularly 
well suited for computations. 

The geometric phase for an evolving quantum system prepared in a pure state 
is undefined if the initial and final states are orthogonal. For quantum systems in 
mixed states a similar situation may occur. In a concluding section we introduce 
higher order geometric phases, and prove that under a generically satisfied condition 
there will always be a well-defined geometric phase of some order. We also relate 
the higher order geometric phases to the off-diagonal geometric phases studied by 
Manini and Pistolesi [26], and Mukunda et al. |27| . 

2. Bundles of purifications 

2.1. Mixed states and purifications. Mixed quantum states can be represented 
by density operators. A density operator is a self-adjoint, nonnegative, trace-class 
operator with unit trace. We write T>(T-L) for the space of density operators on a 
Hilbert space %, and T> n (H) for the subspace of density operators with finite rank 
at most n. 

A state is called pure if its density operator has rank 1. In quantum mechanics, 
especially quantum information theory, purification refers to the fact that every 
density operator can be thought of as representing the reduced state of some pure 
state. Indeed, one can show that if p is a density operator acting on H, there exists 
a Hilbert space K, and a density operator R representing a pure state in H ® K, 
such that p is the partial trace of R with respect to /C. In this paper we only 
consider quantum systems whose states are represented by density operators with 
finite rank, and by the standard purification of the density operators in T> n {'H) we 
mean the surjective map 

(1) ir:S(U®C n *)^V n (H), |*) Tr c „. |*)(*|. 

Here C™* is the space of linear functionals on C", and S(H<E)C n *) is the unit sphere 
in^OC™*. 

Evolving mixed states will be represented by piecewise smooth curves of density 
operators that, for convenience, are assumed to be defined on an unspecified interval 
< t < t. Operators on C will be represented by matrices with respect to the 
canonical basis, whose fc th member we denote by e^. We write 1„ and 0„ for the 
n x n identity matrix and zero matrix, respectively. 

2.2. Bundles of purifications over orbits of isospectral density operators. 

When you purify a state you artificially add extraneous, perhaps inaccessible, in- 
formation to your state. One way to get rid of this extra information is to use 
symplectic reduction. Montgomery [25| showed that if one does symplectic reduc- 
tion on % (g) C n * by the action of the unitary group of C"* , then the reduced phase 
spaces can be identified with the (co) adjoint orbits of isospectral density operators 
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in T> n (H), and the canonical reduced space submersions are given by certain spec- 
tral weighted eigenframe bundles. (A precise reformulation of Montgomery's result 
will be given below.) We show that in this context, there is a natural definition of 
geometric phase for unitarily evolving mixed states which generalizes the definition 
of Sjoqvist et al. But before we do that we need to translate Montgomery's result 
into a language appropriate for our purposes. We begin by identifying T-L<SiC n * with 
the space of linear maps from C" to W, equipped with the Hilbert-Schmidt Her- 
mitian inner product. Then the standard purification ([I]) is given by n(^) — vpv]/t 
A density operator whose evolution is governed by a von Neumann equation stays 
in a single orbit of the left conjugation action of the unitary group U(H) on T>(H). 
These orbits are in one-to-one correspondence with the set of possible spectra for 
density operators on %, where by the spectrum of a rank n density operator we 
mean the nonincreasing sequence a = {p\,P2, ■ ■ ■ ,Pn) of its, not necessarily distinct, 
positive eigenvalues. Throughout we assume a to be fixed, and we write 2?(<x) for 
the corresponding orbit of density operators. Also, we refer to the integer n as the 
length of a. 

Let -P(cx) be the diagonal n x n matrix with diagonal a, and define ir : S(a) — > 
X>(cr) as the restriction of to S(a) = e S(H®C n *) : = P(a)}. Then tt is 
a principal fiber bundle with gauge group U(cr) — {U € U(n) : UP (a) = P(<j)U}. 
We equip S(a) with the metric given by 2h times the real part of the Hilbert- 
Schmidt Hermitian inner product: G(X, Y) = hTr(X^Y + Y^X). 

The vertical and horizontal bundles of 7r are the subbundles V<S(er) and H6>(er), 
respectively, of the tangent bundle TS(a). The former consists of the kernels of the 
differential of 7T, and the latter of the kernels' orthogonal complements with respect 
to G. Vectors in V<S(er) and HS(a) are called vertical and horizontal, respectively, 
and a curve in S(o~) is called horizontal if all of its velocity vectors are horizontal. 
Recall that for every curve p in T>(o~) and every in the fiber over p(0) there is 
a unique horizontal curve 'J in S(a) that starts at 'J'o and projects onto p. This 
curve is the horizontal lift of p extending from ty . Furthermore, if p happens to 
be a closed curve, then the holonomy of p with base point is the unique matrix 
Hol(p, ^o) m U{a) satisfying ^(r) = Hol(p, V^o)- We define the geometric phase 
of p to be the argument of the spectral weighted trace of Ho\(p, ^q): 

(2) 7geo(p) = argTr (P(a) Hol(p, ¥„)) ■ 

Observe that ^ does not depend on because the holonomy of p associated 
with a different base point is conjugate equivalent to Ho\(p, \&o) m U{a), and the 
matrices in U(a) commute with P(cr). 

3. Geometric phase 

3.1. Geometric phase for unitarily evolving mixed states. Inspired by the 
classical work of Pancharatnam [28J, Samuel and Bhandari [29 generalized Berry's 
phase to a relative geometric phase for non-cyclic unitary evolutions of pure quan- 
tum states. Here we extend Q to a relative geometric phase for non-cyclic unitary 
evolutions of mixed quantum states. 

Let p be a curve in T>(o~). The horizontal lifts of p define a parallel transport 
operator H[p] from the fiber over p(0) onto the fiber over p(r) according to nfpj^o = 
^11 (t), where is the horizontal lift of p extending from vPq- We define the 
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geometric phase factor and geometric phase of p to be 

(3) Cgeo(p) = Tr (# o n[p]* ) , 

(4) 7goo(p) = argCgeo(p), 



respectively. Observe that the geometric phase factor does not depend on ^ since 
the right hand side of ^ is invariant under the transitive action of the gauge group 
on the fiber over p(0). Definition Q reduces to Q if the curve is closed. However, 
for this definition to be useful in practice, we need a way to derive a horizontal lift 
from an arbitrary lift of a curve in T>(o~). For this we need a connection form. 

The Lie algebra of the gauge group is u(ct) = {£ G u(n) : £P(cr) — P and 
the infinitesimal generators of the gauge group yield canonical isomorphisms be- 
tween u(fi) and the fibers in VS(a): u(a) 3 £ i-> ^ G V*<S(er). Further- 
more, the horizontal bundle is the kernel bundle of the mechanical connection form 
A : TS(er) -> u(cr) defined by = I^J*, where I : S(<r) x u(a) -> u(cr)* and 
jT" : T5(cr) — > u(cr)* are the locked inertia tensor and moment map, respectively: 

If /; is a curve in T>(a), and is any lift of p to <S(cr), then the horizontal lift of p 
that extends from 11/(0) is 

(5) ¥||(t) =V(t)exp+(-J Av(*)dt\, 

where exp + is the positive time-ordered exponential. The exponential factor in ([5| 
has a nice geometric interpretation; it describes a curve in the gauge group 14(a) 
whose continuous action on S(a) pushes ^ onto a horizontal curve, see Figure [I] 




FIGURE 1 . Action of the exponential factor in Equation ^ . 



Equations Q and (|5| give a truly general definition of geometric phase for 
unitary evolutions of mixed states. We will now derive an explicit formula for 
the geometric phase, suitable for computations. Thus, let mi, m.2, ■ ■ ■ mi be the 
multiplicities of the distinct eigenvalues in a, where mi is the multiplicity of the 
greatest eigenvalue, m,2 of the second greatest eigenvalue and so forth, and define I 
diagonal matrices by 

P^k diag (0 mi , . . . , mfc _ 1 , l mfc , Orrik + i 3 • • • ) ^m; ) • 
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A straightforward computation shows that I*Efc XEf-P{cr) 1 )-n = J'^(X)-r) 
for every vector X tangent to S(a) at if? and i] in u(er). Hence, 

(6) A^{X)=Y J E k ^XE k P{a)-\ 

k 

Moreover, if U is a time-evolution operator that generates p, i.e., p(t) — U (t) p(0)U (ty , 
we can take ^ = U^q, where ^0 is an y element in the fiber over p(0). According 
to {§,([5]), and {6}, 

lgco {p) =argTr U^r)^ exp + ^- jf £ ^^(iM^o^/V)' 1 At 



k / / 

3.2. The geometric phase of Sjoqvist et al. Generically, a density operator of 
a mixed state is nondegenerate. In [20], the authors gave a definition of geometric 
phase for parallel transported nondegenerate mixed states. We will quickly review 
it here, and then show that it is a special case of our definition. 

Let the curve p = X^-Pfcl^X^I represent a unitary evolution of a nondegenerate 
density operator. Following Sjoqvist et al. we say that the initial density operator 
p(0) is parallel transported provided all of its eigenkets are parallel transported: 
(k\k) = 0. The geometric phase of p in the sense of [20 is then the argument of the 
expected value of the relative phase factors for the eigenket curves: 

(7) 7g co(p)=arg(^p fc (fc(0)|fc(r))). 

fc 

To show that Q generalizes ([7]), let p be any curve in T>(a), Vf^ be any horizontal 
lift of p, and set Wo = ^||(0)- Define curves of normalized, and at each instant 
pairwise orthogonal, kets |fc) = p k ^ffy^ek- Then p = ^2 k Pk\k)(k\, and it follows 
from the horizontality of W||, and the observation (k\k) = e\Aq,^ (W^e/., that each 
curve of eigenkets is parallel. (In fact, using one can show that the horizontality 
condition Ay., (^\\) = is equivalent to the statement that, at each instant t, the 

velocity vector \k(t)) is orthogonal to the eigenspace of p{t) corresponding to the 
eigenvalue p k .) Now 

Tr (*Jn[p]* ) =J2 e i^\\m\\(r)e k = £p fc <fc(0)|fc(r)>. 

k k 

This shows that our definition of geometric phase generalizes the definition of geo- 
metric phase in |20j . 

3.3. A mixture of neutrons in a magnetic field I. Consider an ensemble of 
neutrons subjected to an external magnetic field, and assume that the ensemble 
is prepared so that the proportion of neutrons with spin up polarization is p\ and 
the proportion with spin down polarization is p 2l where p\ ^ p 2 - The initial state 
of the spin part of the ensemble's wave function can be represented by the density 
operator 

Pi 0" 

P2 ■ 

(Here we have modeled the system on C 2 in such a way that ei and e 2 represent the 
spin up and spin down states, respectively.) If we further assume that the magnetic 
field is static and uniform with direction n = (sin 6, 0, cos 8) T , then the Hamiltonian 
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G 
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is given by H = -un • cr, where u> depends on the magnitudes of the magnetic field 
and of the neutrons' magnetic moments, and cr is the vector of Pauli matrices. The 
time-evolution operator associated with H is the one-parameter family of unitarics 



U(t) = cos(wi) 



1 
1 



i sin(wt) 



cos ( 
sin t 



sin 6 
- cos 6 



Observe that this family is non-parallel transporting if cos 8 7^ 0. Now, a lift of the 
ensemble's evolution curve p = Upoll 1 ' is given by VP = U^o, where 



*n = 



Pi 






P2 



If the final time equals n/u>, the curve p returns to its initial state, and the geometric 
phase of p is 

7g eo0) = argTr(M'JC/(7r/a;)*oexp + f- / ' A*(i>)dt) 



■(■ 



arg p x e 



Z7r(l+COS t 



P2e 



,7r(l — c 



4. Higher order geometric phases 

If the initial and final purifications and n[p]\E r o are orthogonal, the geometric 
phase factor Q vanishes, and, consequently, the geometric phase Q is not defined. 
In this section we show how submatrices of \&qII[p]^o can be combined into higher 
order geometric phase factors that might be nonzero even if the geometric phase 
factor is zero. For a cyclically evolving state these phase factors will be formal 
homogeneous polynomials in the matrix entries of the holonomy of the evolution 
curve (hence the name 'higher order' geometric phase factors). We also formulate 
a fairly general condition that guarantees the existence of a well-defined geometric 
phase of some order not greater than the length of cr, and we observe that this 
condition is satisfied by cyclic evolutions. 

4.1. Geometric phase of degree d. Suppose p is a curve in T>(a), and is an 
element in the fiber over p(0). For each sequence J = (ji, j2, ■ ■ ■ ,jd) of elements 
from {1, 2, . . . , /} we define the geometric phase factor and geometric phase for p 
corresponding to J by 
d 

UP) = Tr n E jk ^ U[p]^ E jh+1 ( Jd+1 = ii), 7j (p) = argO(p), 



k=l 



respectively. The geometric phase factor (J3j) can be recovered from the geometric 
phase factors corresponding to single element sequences: Cgeo(p) = Y^j=iC(j)(p)- 
More generally, if d is any positive integer, then Tr((\E'on[ / o] v I'o) ci ) = SjO(/°)j 
where the sum is over all possible sequences J = (ji, j2, ■ ■ ■ ,jd) that can be formed 
from the elements in {1, 2, . . . , I}. We define the d th order geometric phase factor 
and geometric phase of p as 

Cgco (P) = ^ O (P) , 7gco (P) = arg Cgoo (p) ■ 
J 

Below we show that there is always a nonzero geometric phase factor of an order 
not greater than the length of a if ^Jnfpj^I'o h as a nonzero eigenvalue. 
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For cyclic evolutions the geometric phase factor (j(p) vanishes if the sequence J 
is not constant. Indeed, 

f[E jh Hu[^ E jh+1 = (WHoKp^o^if h =h = -=U, 
jJJ 3k 3k+1 \0 n otherwise, 

due to the block diagonal structure of the matrices in U(<j). It is apparent from 
this equation that (j(p) is a formal homogeneous polynomial of degree d in the 
entries of Hol(p,* )- Moreover, the matrix *JlI[p]vl/ = P{a) Hol(p, has a 
nonzero eigenvalue because it represents an isomorphism of C™. Thus, a cyclic 
p has a well-defined geometric phase of some order. (On the other hand, if p is 
open and the initial and final states p(0) and p{j) are distinguishable [301 13 lj . then 
\E , gII[/>]\I f o = 0„ why (j(p) = for every sequence J, c.f. |32].) 

4.2. A mixture of neutrons in a magnetic field II. Consider the ensemble of 
neutrons in section |3.3[ but now assume that H is such that 

U(t) 



cos(wt) sin(w£) 
sin(W) — cos(wi) 

If the final time is it/2uj, then Cgeo(p) = and C|eo(/ 5 ) = 2??iP2- Consequently, the 
geometric phase of p is not defined, but its second order geometric phase is defined. 

4.3. The existence of higher order geometric phases. If ^ , JiI[p]\E'o has at 
least one nonzero eigenvalue, then p has a well-defined geometric phase of some 
order less than or equal to the length of a. To see this observe that the characteristic 
polynomial of \I r QlI[p] 1 I r o can be written 

n 

det (A - *Jn[p]* ) = ^(-l) fe Sfe A"- fe , 

k=0 

where 



3 = 1 

Therefore Cg 0O (p) + for somc d < n if det (A - *J,n[p]* ) = for some nonzero 
A. 

4.4. Relations to the off-diagonal geometric phase factors of Manini and 
Pistolesi. It is customary to represent an evolving density operator as an incoher- 
ent sum of 1-dimcnsional time-dependent projectors defined by the normalized kets 
in a comoving eigenframe: p(t) = Pk\k(t)) (k(t)\. This representation is unique 
for nondegenerate mixed states, and if the eigenkets are parallel transported, the 
geometric phase factor is the expected value of the relative phase factors (fc(0)|fc(r)), 
see section [3721 

If all the relative phase factors vanish, the geometric phase is not defined. Manini 
and Pistolesi |26| then proposed to study the gauge invariant off-diagonal phase 
factors 

7g fe3 ,... )fc » = n<Mo)i^ + i(T)> G**fi) - 

Mukunda et al. [57] extended the work of Manini and Pistolesi. They expressed 
the off-diagonal phase factors in terms of Bargmann invariants, and used the fact 
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that the Bargmann invariants satisfy certain relations, to give a full description of 
how the off-diagonal phase factors are interrelated. 

For degenerate mixed states the off-diagonal phase factors need not be gauge 
invariant. However, appropriate combinations of them are. For example, 



where each index k s run through all integers satisfying < k s — Y^Li nii < frij s . 
Hence, using Q and the results of Mukunda et al., one can, in principle, describe 
all relations between the higher order geometric phases. 



We have presented a fiber bundle framework in which geometric phase for mixed 
quantum states has a natural definition. This definition extends the standard defi- 
nition of geometric phase of Sjoqvist et al. |20| . which is based on interferometry, 
and is particularly well suited for computations in quantum physics. Furthermore, 
we have shown that the geometric phase is the first in a sequence of higher order 
geometric phases, and we have formulated a criterion that guarantees the existence 
of a geometric phase of some order. This criterion is met by cyclic evolutions of 
mixed states. 

Another approach to geometric phase for mixed states is given in Uhlmann [18 
IT!?] . Uhlmann's approach only applies to mixed states represented by invertible 
density operators, and according to Slater [33 it generally yields different outcomes 
than that of Sjoqvist et al., and hence ours. The discrepancy is probably due to a 
difference in the notions of parallel transport [34], [32] . Also Singh et al. [351 give a 
definition of geometric phase for parallel transported mixed state^] and the paper 
by Chaturvedi et al. |36| contains a differential geometrical approach to geometric 
phase for mixed states. Comparative studies of our definition and those in [35J and 
[36 will be published in forthcoming papers. 

Furthermore, in a forthcoming paper, off-diagonal geometric phases for mixed 
states will be defined and investigated. Such phases have been studied by Filipp and 
Sjoqvist [371 1351 133] . In [33], the authors defined "off-diagonal quantum holonomy 
invariants" for the bundle of purifications by Uhlmann, and compared these with the 
off-diagonal geometric phase factors defined in [38 . The conclusion was that there is 
"a general discrepancy for these two approaches related to a fundamental difference 
in the treatment of parallel transport of quantum states". The off-diagonal quantum 
holonomy invariants can be similarly defined for the bundle of purifications in the 
current paper, and we will investigate how these are related to the off-diagonal 
geometric phases by Filipp and Sjoqvist. 



1 Unfortunately, the proof of gauge invariance of the proposed geometric phase in 1351 lacks 
vital parts. The authors are forced to introduce certain initial values in the proof that has as 
a consequence that the proof does not imply full gauge invariance. Moreover, the definition of 
geometric phase is such that for an evolution of a noninvertible density operator, the holonomy 
of an artificially chosen frame of kets in the eigenspace corresponding to the density operator's 
zero eigenvalue seems to contribute to the geometric phase. We don't think this is correct from a 
physical point of view. Also, the definition does not, in an obvious way, reduce to the Aharonov- 
Anandan geometric phase for pure states, and it is not clear how to generalize the definition to 
infinite dimensional Hilbert spaces. 



(8) 
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